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Abstract. We study the nonlinear inhomogcncous wave equation in one space dimension: 
vtt — T(v, x) X x = 0. By constructing some "decoupled" Riccati type equations for smooth solutions, 
we provide a singularity formation result without restrictions on the total variation of unknown, which 
generalize earlier singularity results of Lax and the first author. These results are applied to several 
one-dimensional hyperbolic models, such as compressible Euler flows with a general pressure law, 
elasticity in an inhomogeneous medium, transverse MHD flow, and compressible flow in a variable 
area duct. 
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1. Introduction. In this paper, we consider the initial value problems for the 
second-order quasilinear inhomogeneous wave equation in one space dimension, 

v tt ~ T(v, x) xx = 0, (1.1) 

where (x,t) € R x R + , v(x,t) G R, and T(v,x) is a smooth function satisfying 



These assumptions imply that the equation is hyperbolic and genuinely nonlinear, 
so that solutions exhibit wave-like behavior. Signals propagate in a forward and 
backward direction with local nonlinear wavespeeds 

c = yTv an d — c — — , 

respectively. Equation (|1.1[) includes a wide variety of interesting physical systems, 
such as one-dimensional nonlinear elasticity in an inhomogeneous medium [5] [6] . In 
the context of elasticity, x is a material coordinate, v is the strain and T is the elastic 
stress. 

It is well known that solutions of nonlinear hyperbolic equations generally form 
shock waves in finite time. Shocks form as a result of gradient blowup, which is a 
consequence of genuine nonlinearity. In this paper, we study smooth solutions, and 
we are particularly interested in the lifetime of such smooth solutions. We study the 
dynamical system which governs the growth of gradients, and use this to give esti- 
mates on the lifetimes of smooth solutions. This was carried out for the homogeneous 
nonlinear wave equation, obtained by taking T = T(v), by Lax in [TJ, and exten- 
sions to larger systems (with restrictions on the data) were obtained by John [5] and 
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It is convenient to write as a system, by setting 



u(x, t) = 



v t dx, p(v,x) = -T, 



which yields the first-order system 



v t -u x = 0, 
u t + p(v,x) x = 0, 



(1.2) 
(1.3) 



with 



p v < 0, p vv > 0. 



(1.4) 



We make the further assumption that p — p{v, x) is a C 3 function, and that bounded- 
ness of derivatives of p and the sign conditions (jl.4l) are uniform in x. In particular, 
p x and p xx are uniformly bounded in x. 

For smooth solutions, equations (|1.2p . (|1.3|) model compressible inviscid flow for 
general pressure laws in material coordinates (the Euler equations in a Lagrangian 
frame). Here p is the pressure, v is the specific volume, and u is the fluid velocity. A 
general pressure is described by p — p(v, S), but the entropy S = S(x) is stationary 
as long as the flow is smooth, by the Second Law of Thermodynamics. Moreover, 
smooth solutions of the 5x5 system of one-dimensional transverse flow in Magne- 
tohydrodynamics (MHD), which models a fluid coupled to a transverse (magnetic) 
vector field, and Eulerian flow in a variable area duct can also be modeled by (|1.2[) , 



For homogeneous isentropic flow, we have p = p(v), and (|1.2[) . (| 1 . 3[) is known as 
the the p-system. This is a simplified 2x2 system which admits a coordinate system 
of Riemann invariants, and for which Lax proved that any nontrivial data will form a 
shock wave in finite time [7]. F. John and T.-P. Liu extended Lax's results to general 
systems of conservation laws of the form 



with u — (ui,--- ,u n ) £ R™, for n ^ 1 and (x,t) £ R x R + . John proved that 
gradients will blow up for small compactly supported data, and this was later gener- 
alized by Liu, provided the total variation of the data is small enough [5j [10] . There 
are also some singularity results for multi-dimensional conservation laws, subject to a 
restrictive "null condition" |T2l [13] . In a recent paper [1] , the first author generalizes 
the singularity formation results in [7] to the 3x3 compressible Euler equations with 
polytropic ideal gas. 

In this paper, we give singularity formation results for (|1.2[) . (|1.3[) . or equivalently 
(ll.ip . In [5][TU], small solutions are expanded asymptotically along integral curves, 
and wave interactions are treated quadratically. Here we take a different point of 
view, considering three wave families, namely forward and backward waves, as in the 
p-system, coupled with stationary waves, which carry entropy variations and inhomo- 
geneity. We are then able to treat all waves without regard to their wave strength. 

Our first task is to define the rarefactive and compressive character (R/C char- 
acter) of the nonlinear (non-stationary) waves in smooth solutions. In a 2 x 2 system, 
which is diagonal when expressed in Riemann invariants, it is clear when a wave is 
rarefying or compressing. However, in larger systems, this distinction is not clear as 



«t + f(u)x = 0, 



(1.5) 
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waves of different families generally interact continuously and cannot be decoupled . 
The R/C character is a quantitative measure of how much rarefaction or compression 
is in the solution at any point. 

In stationary solutions (including those with contact discontinuities), which in- 
clude no compression or rarefaction, the pressure p and velocity u are constant, by 
the Rankine-Hugoniot conditions. We therefore use changes in p (or u) to define the 
R/C character of the solution. In doing so, we take the directional derivative along 
the opposite characteristic, to minimize the effect of waves of the opposite family, see 
PP. We use the superscripts "\" and "/" to denote the directional derivatives along 
backward and forward characteristics, respectively, so that 

y = d t -cd x and ' = d t + cd x , (1.6) 

where c = \/—p v is the (local) wavespeed in Lagrangian coordinates, c.f. [7]. 

Definition 1. // the solutions of il.ty) . il.3\) are smooth in an open set U of 
the (x,t)-plane and A is a point in U , then we say the solution is forward (backward) 
rare] "active at A, if and only if p* < (p 1 < 0); it is forward (backward) compressive 
at A, if and only if p* > (p' > 0). 

It is convenient to introduce the change of variables 

r v * 

h(v,x)= / \J — p v dv, and [i = x, (1-7) 



where v* is a constant or infinity. Then by calculating p y and p' , we introduce equiv- 
alent variables a and /3, defined by 

a = — 7 = u x +h x + , 
c z c 

a- P' u Pv- 
p = — k =u x - h x . 

cr c 

Here u ± h are the Riemann invariants for the corresponding isentropic p-system. 
Thus a and /3 are direct generalizations of the derivatives of the Riemann invariants. 
For smooth solutions, we derive Riccati type ODEs for a and j3, which provide a 
framework for studying smooth solutions and gradient blowup. 
Theorem 1. In (Op , smooth solutions satisfy 

a' = -^) h (3a + (3) + ^(al3-a 2 ), (1.8) 
2 c 2 

P = c ^) h { a + m + c ^{*f3-p 2 ), (i.9) 

2 c 2 

with Ch > 0. 

Equations (jl.8l) and (|1.9[) can be decoupled by use of an integrating factor. Define 

y=^/ca-I = ^/c(u + h) x + ^-I, and (1.10) 

q = ^ + I = ^( u -h) x -^+I, (1.11) 

\ c 



where 



I = I(h,n)= t^c^udh, (1.12) 

An 2 C 
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and h is a constant. 

Theorem 2. For smooth solutions of M.2\) , U.S\) , we have 



where 



y = a a + a 1 y-a 2 y 2 , (1-13) 
q = a - ai q - a 2 g 2 , (1-14) 



ao = -d, + lVc(^)hP, ~ c(^) h I - ^I 2 , (1.15) 

Z C C Z*\J c 



ai = -(2y/il) h , (1.16) 
a 2 . J* > 0. (1.17) 



We note that these are not closed ODEs, since both the directional derivatives and 
the coefficients are dependent on the underlying solution of (|1.2p . (|1.3|) . Nevertheless, 
as in [7J, we are able to compare them to closed ODEs and derive bounds on the 
lifespan of smooth solutions. For convenience, we only consider smooth, i.e. C°°, 
initial data, although our results also apply to C 2 initial data. 

Theorem 3. Assume the smooth solution takes values in some compact set K,, 
uniformly in x. Then there is some constant N > depending only on K,, such that, 
if y or q is less than —N somewhere in the initial data, then \u x \ and/or \v x \ blow up 
in finite time: that is, there is some critical time T* < oo, so that 

max { | u x | , | v x | } — » oo as t — > T» . 



Blowup of the gradient does not mean that the solution fails to exist: rather, 
this usually heralds the formation of a shock, and the associated decay of solutions. 
Generally, solutions are continued as weak solutions which contain shocks. However, 
in order to study shock propagation our system needs to be in conservation form, and 
we do not consider those issues here. 

Our assumption that the solution stay in a compact set K, can be stated technically 

as 

\h\ <Bi, A 2 <c< fl a , A 3 < c h < B 3 , (1.18) 
|c M | < B A , |c MM | < B 5 , \c hfl \ < B 6 , |p M | < B 7 , \pnp\ < B s , 

where Ai and Bi are all positive constants. These conditions simply mean that we do 
not leave the domain K, in the solutions we consider, and allow us to stay away from 
vacuum and other points where the equation may become singular. From a practical 
point of view, the only restriction on the data (apart from smoothness conditions) is 
that the solution stays away from vacuum: in general we do not expect the pressure p 
to grow unboundedly. For the p-system, an invariant region provides an upper bound 
for the pressure, and although there is not an invariant region for the full system, 
physically we expect the pressure to remain finite. We note that the assumption that 
the solution does not form a vacuum is implicit in jjj, but this is reasonable as a 
vacuum cannot form in finite time, see |17[ 118]. 

Theorem [3] implies that gradients of solutions blow up if the initial compressions 
are strong enough. When the variation of entropy is mild, N is close to zero, so the 
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shock free solutions are "almost rarefactive" , which is consistent with Lax's singularity 
formation results in [7]. In [Tl] and a forthcoming paper [2J, examples of solutions 
containing compressive waves are constructed, but the gradients of those solutions 
remain finite. 

In a series of recent papers [14j [15l [16] , the possibility of time-periodic solutions 
in the compressible Euler equations has been demonstrated. A critical feature of this 
study is how the R/C structure of waves can change across a contact discontinuity. 
In this paper we check the consistency of our results, which presuppose a smooth 
entropy field, with those results. By studying the R/C structure and the way in 
which it can change further, we expect eventually to see time-periodic solutions as in 
[T4l H3] with piecewise smooth entropy, consisting of both contact discontinuities and 
smooth entropy variations. 

Our results apply directly to a number of systems having structure similar to 
the inhomogeneous nonlinear wave equation. In addition to equations (|1.2p . (|1.3I) . 
we apply these ideas to transverse flow in one-dimensional magnetohydrodynamics 
(MHD), a 5 x 5 system modeling fluid motion coupled to a transverse magnetic field. 
By restricting to a polytropic ideal gas, we are able to give a stronger singularity 
formation result than that of Theorem 02 similar to that of [1] . Finally, we consider 
inviscid compressible flow in a varying area duct, to which similar ideas apply. 

The paper is arranged as follows. In Section 2, we give the background of the 
equations and establish some useful identities. In Section 3, we define rarefactive and 
compressive waves. In Sections 4 and 5, we prove Theorems HH21 respectively. In 
Section 6, we consider the consistency of R/C structures, and in Sections 7 and 8, we 
apply our ideas to one-dimensional transverse flow of MHD and compressible flow in 
a varying duct. 

2. Coordinates and background. We focus on equations (jl.2l) . (|1.3p from 
now on. We assume that p(v, x) is a smooth function of v and x, satisfying 

p v < 0, and p vv > 0, 

for all x and v = v(x,t) £ (0, oo). These conditions imply hyperbolicity and genuine 
nonlinearity, respectively. The vacuum state corresponds to v = oo; however, since we 
assume the data remains in a compact set, we will not address questions at vacuum. 
The local absolute wavespeed is 

c(v,x) = v^>0. (2.1) 
We make the change of variables 

h(v,x)= / c dv = \/—Pv dv, 

J v J V 

see (|1.7j) . where v* > is a convenient constant (or oo if the integral converges 
uniformly). Since p(v,x) is smooth, the function h(v,x) is also smooth with respect 
to v and x. Moreover, since 

h v (v,x) = -c< 0, (2.2) 
the inverse function v = v(h, /j) is smooth with respect to h and /z, where we have set 



/! = x. 



(2.3) 
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For any function f(v, x), we will write 

f(h,n) = f(h(v,x),n) = f(v,x), 
without ambiguity, and we use the subscript notation 

f x = £f(v( x ,t),x), f T =£f(v,x), f v = JLf(v,x), 



for the various partial derivatives of /. 

We can relate the different partial derivatives as follows: by 

Vh = -~, so that v t = v h h t = --h t . (2.4) 

c c 

Furthermore, since 

- c 2 = p v = p h h v = -cp h , we have p h = c, (2.5) 

and so 

Px = Vhhx +P^ = ch x + Pp. 
Thus, for smooth solutions, (|1.2j) . (|1.3|) can be written as 



h t + cu x = 0, (2.6) 

u t + ch x +p M = 0. (2.7) 

Next, differentiating v = v(h(v, x), /i) with respect to x, and recalling (|2.3|) and 
(f2T4|) . we get 

= Vhh— + v^, so that = —Vhh— — —. 
Similarly, for any function f(v,x), we get 

= fvVfi^r fx 

= -h-+h, (2.8) 
c 

and 

fh = fvV h = --• (2.9) 
c 

It follows that, for smooth solutions, p and c are smooth with respect to h and /i. 

Moreover, using (|2.8[) and (|2.9j) . any quantities such as the R/C character, ODEs, 
and singularity formation results in this paper can all be expressed in the variables 
(v,x) instead of (h,fx). 
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3. Compressive and Rarefactive waves. In this section, we define the rar- 
efactive (R) and compressive (C) characters of (ll.2[) . (I1.3[) . which quantitively indicate 
the amount of rarefaction or compression in the solutions at any point. 

We first consider isentropic flow, for which p = p(v) and our system reduces to 
the p-system, 

h t + cu x — 0, (3-1) 
u t + ch x = 0. (3.2) 

The Riemann invariants s — u + h and r = u — h satisfy the diagonal system 

s t + cs x = 0, r t - cr x = 0, 

and so are constant along forward and backward characteristics, 

dx dx 

dt = +C ' dt = ~ C ' 

respectively. In an isentropic domain, because the system is diagonal, it is clear 
when waves are compressive or rarefactive: indeed, the amount of compression or 
rarefaction can be measured by derivatives of the appropriate Riemann invariant. 
There are several equivalent conditions, and for us it is convenient to consider the 
change in pressure: if p decreases as we traverse the wave from front to back, the 
wave is rarefactive (R), while if p increases, it is compressive (C) [TJ [71 [17] ■ This 
is consistent with the entropy condition for shocks, which states that the pressure is 
always larger behind a shock. 

When p = p(v,x) explicitly depends on x, we first consider stationary solutions, 
in which there are no compressive or rarefactive waves, so the R/C characters should 
vanish. In stationary solutions, the pressure p is constant, 

Pt = Pv v t = and p x = -u t = 0, 

so its directional directives are zero. Physically, in gas dynamics, this means that 
pressure is not impacted by the variation of entropy, c.f. [T]. Thus, by considering the 
directional derivatives of p along the opposite characteristics, we obtain Definition [T] 
the solution is forward (backward) rarefactive at A, if and only if p y < (p' < 0); it 
is forward (backward) compressive at A, if and only if p~ > (p' > 0). 

We could also use u to define the R/C character, as this is also constant in 
stationary solutions. We define 

a = — r and p = — -. 

c z c 

Lemma 1. For smooth solutions of hl.2\) , H.3\) , we have 

p' = —cu and p^ — cu , 

while also 

a = u x + h x + — and (3 = u x — h x — —. (3-3) 
c c 
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Proof. By (O, (T3) and (ETD) . 

cu' = cu t + c 2 u x = -cp x - p v v t = -p', 

and similarly p x = cu' . By (I2.5[) and (I2.6[) , 

-c 2 /3 = p' = p t + cp x 

= p h h t + c(phh x +p f _ l ) 

= -c 2 {u x - h x - — ), 
c 

and similarly for a, so (13.31) follows. □ 

Corollary 1. The R/C character of a smooth solution is given by: 

Forward R iff a > 0, 

Forward C iff a < 0, 

Backward R iff ft > 0, 

Backward C iff [3 < 0. 

Moreover, provided the solution values remain in K,, 

\a\ or \[3\ oo iff \u x \ or \v x \ -> oo. (3.4) 



Proof. Clearly, by (l3~3| . 

p x ^ <S> a ^ 0, and p' ^ [3 ^ 0. 

By Lemma [T] 

a + f3 = 2u x , a-[3 = 2(h x + ^), 

c 

and (13. 4[) follows since p M and c remain finite. □ 

In an isentropic domain, i.e. p = p(w), it is clear that 

s x — a and r x = /3, 

so we can regard a and /3 as direct generalizations of the derivatives of the Riemann 
invariants. 

4. Differential equations for gradients. In this section, we consider the char- 
acteristic decompositions of smooth solutions. By considering the directional deriva- 
tives of a and /3, we derive the ODEs for a and f3 as stated in Theorem [T] 

Proof of Theorem^ We show (|1.8[) . since (|1.9I) follows in exactly the same way. 
We have 

a = (u x + h x + ^-) t + c(u x + h x + —) x 
c c 

= (u xt + ch xx ) + (h xt + cu xx ) + ( — ) t -I- c( — ) x (4.1) 

c c 

= (u t + ch x ) x + (h t + cu x ) x - c x h x - c x u x + ( — ) t + c(-^) x . 

c c 

By (nu, (mm, 
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(ht + cu x ) x = 0. 
Thus the right hand side of (|4.1[) is 

-Pfihhx -p^ - {c h h x + c fl )(h x +u x ) + (-^) h (h t +ch x ) + c(-^) M , 

since f x = f h h x + f u for any function /. 
By ((33)) and (|2T6l) . we have 



h x = a — u x — — , and h t = —cu x , 
c 



and by (l23j) . 



c( — )/j = c„ - c ft — . 

c c 



Thus (|4.2[) can be simplified to 

-P^h {a-u x - -^) - pan - {c h a - c h u x - c h — + c„)(a - — ) 
c c c 

+ ( — ) h (ca - 2cu x - p u ) + c(— )a 
c c 

= -c(-^) h u x + [c h u x - c(—) h ]a - c h a 2 . 
c c 

Finally, making the substitution 

a + (3 

M.T = ~ 



yields (TL81) . 

By ([LID, (H3J and (l23j) . 

> W-Pv)v = c v = c h h v = Ch(-c), 
so > 0, and the theorem is proved. □ 

Corollary 2. For smooth solutions of U.ty) , U.3\) , we have 

«' = -j(-),(3a + /?)-^(a/3-« 2 ), 
4 p v 2c 

and 

P = j(-Ua + 3/3) -^(a/3-/3 2 ). 
Proo/. By flZIJ), (EH), ([21]) and (US), we have 
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and by fO}, 



and the corollary follows from Theorem [TJ □ 

We now make another change of variables and transform (|1.8I) and (|1.90 into 
decoupled differential equations by use of an integrating factor. 

Proof of Theorem^ First, the condition ai > follows immediately from (|4.4[) . 

By and flUD, 

T) 

h' = h t + ch x = —cu x + ch x = —c(/3 H — — ), 



so that 



Hence, (I1.8P can be written 



ti 



_ _ Ea 

c c 



, C Pf,, ti PfJ , C h ti p M C h 2 

a = --{ — )h(3a £) + —a( £) - — of. 

2 c c c 2 c c 2 

We move the terms including ti to the left hand side, so 

a - o(—)hh + —ah 
2 c 2c 



Now, by (|2"3j) . 



1(^ + ( _^ K _£,,„)«, (4.6) 



- - _„„ = -c(-) h - — = -c(-) h - —, 



and, since p! — c, we have 



Vca + -^—j=ah + —c^a = (y/ca) . 
2\l c 2 



Thus, multiplying (14.61) by y/c and simplifying, we get 
2 c 



\M P -^)hP, - cV~c(^) h a - ^a 2 . (4.7) 
2 c c 2 



y + I 



a 



where / is defined in (|1.12p and satisfies 

I' = I h ti + cI, L 



SINGULARITY FORMATION FOR THE WAVE EQUATION 



11 



Using these in (I4.7[) . we get 

?/ = a + a x y - a 2 y 2 , (4.8) 

where Oq, ai and a 2 are defined in (|1.15l) - (|1.17p . 

The derivation of the differential equation along backward characteristics, 

q x = a - aiq - a 2 q 2 , (4.9) 

where q is defined in (|l.lll) . is similar, and the proof is complete. □ 

Corollary 3. A singularity (gradient blowup) forms if and only if 

\y\ or \q\ oo iff \u x \ or \v x \ oo, 

provided the solution takes values in the compact set K,. 
Proof. By (|TTTTjj) and (1TTT1) . 

y + q= Vc(a + (3), 
y-q= ^~c{a - f3) - 21. 

By (|4.3|) . (| 1 . 1 2[) and compactness, / remains finite, and the result follows. □ 

Because the coefficients do, a%, and a 2 in Theorem [2] don't include derivative 
terms v x , u x , v t or u t , they are lower order when compared to y and q. Using (|2.8j) 
and 1|2.9|) , these coefficients can be expressed in terms of v and x rather than h and 
H, as in Corollary [5] 

In the p-system, p = p(v) and p^ — 0, so that ao = a\ = 0. In this case, (|1.13|) 
and (j!.14|) become 

y = -a 2 y 2 , q' = -a 2 q 2 , (4.10) 

which are exactly the ODEs derived in [7j [8] for the isentropic homogeneous case. 

5. Formation of singularity. We now consider the formation of singularities, 
which take the form of the blowup of gradients u x and/or v x , and correspond to shock 
formation in conservative systems. 

We will study the equations (|1.13[) . (|1.14[) as a dynamical system, even though 
this is not a pure system of ODEs. First, consider the ODE 

where tp± is defined by 

Tjj±(0 = a ± ai t-a 2 e, (5.1) 

with a 2 > 0, and the are treated as constants. The equilibria, if they exist, are the 
roots of the quadratic equation 

a Q ± ai£ - a 2 £ 2 = 0, 

and we have £ > between the roots, and f < otherwise. A typical phase line is 
shown in Figure UTTl 

Whenever there are real roots, say £i < £2, the region {£ > £1} is an invariant 
region for the ODE. In particular, if ap > 0, then the roots are of opposite signs and 
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Fig. 5.1. Phase line fori, = V>±(£). 



the region {£ > 0} is invariant. Moreover, the region {£ < £1} is also invariant, and 
solutions that originate in this interval have only a finite time of existence: it is this 
that drives the gowth and blowup of gradients in the full system (|4.8p . (14.91) . 

For the p-system, we have = 0, so ao = a% = and 1 = 0. In this case, y = 
\fcs x and q = \fcr x are (multiples of) the gradients of the Riemann invariants. By the 
above discussion, the regions {y > 0}, {q > 0} are invariant domains for the system 
(|4.10p . as are the regions {y < 0}, {q < 0}. Thus if y or q is negative somewhere, 
the negative quadratic functions in (|4.10p drive blowup of y or q in finite time, as 
shown by Lax in [7]. From our point of view this is a trivial case, for which we have 
uniform estimates. In general, the forward and backward waves interact nonlinearly 
with the varying stationary background, and we do not expect uniformity. In fact, 
interactions can cause waves to change their R/C character, as demonstrated in [14] 
and the upcoming paper [2]. We expect that a complete analysis of the dynamics of 
(j!.13[) . (|1.14l) will yield a rich variety of new and unexpected phenomena. 

We now prove the breakdown results of Theorem [3] by studying (|4.8j) , (|4.9j) as a 
dynamical system. Our aim is to describe sufficient conditions which imply that the 
gradient blows up in finite time. 

Proof of Theorem [3J Fix a constant 0<v<l, and define 

V4(0 = ao±ai£-(l-^a 2 £ 2 , (5.2) 

so that our ODEs ([478]) . ([479]) can be written 

y' = 4> u + {y) -va 2 y 2 and q s = ijj v _{q) - v a 2 q 2 '■ (5.3) 

Now let N = N(v) < be a uniform lower bound for the (real) roots of i/)±, or 
N = if there are no real roots. Then since a 2 > 0, we have 



V4(£)<0 for every £ < N. 



Now suppose there is some xq such that the data satisfies 

y = y(x o ,0) < N. 
Then for the forward characteristic emanating from (cco,0), we have 

y' < -va 2 y 2 , 

so that the solution satisfies 

1 1 f* 

> — - + v \ a 2 dt, (5.4) 



y(t) - 2/(0) Jo 
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where the integral is taken along the forward characteristic. Since y(0) < and a-i 
is uniformly positive, there is some finite T* such that the right hand side of (|5.4[) 
vanishes, and so we obtain 

y{t) — > — oo as t — > 

A similar calculation holds for q. 

Finally, we calculate the lower bound N. The roots of ip± solve the quadratic 
equations 

r±(O=a ± ai ^-(l-u)a 2 e =0 

and so are 

±ai ± -s/A 2 
£ = — — , where A = a 1 + 4 (1 - v) a a 2 . 

2(1-1/) 02 

The minimum of these is clearly 

|aij_+vg 
2(l-i/)a 2 ' 

and the lower bound N is obtained by maximizing the ratios |ai|/a2 and ao/a2, so 
we look for upper bounds for ao and |ai|, and a lower bound for a-i. 

From the expressions (|1.15[) - (|1.17|) and (|1 . 12|) . we obtain the bounds provided c 
and Ch are bounded away from zero, while the quantities 

h, c, c h , p^, p^ h , Pnhh, and p^h, 

remain finite, and recalling (|2.5p . the bound follows. 

Using ()2.8|) and ()2.9|) . the bounds can also be expressed by the bounds on the 
derivatives of p (and c) with respect to v and x. □ 

When ao > 0, the equilibria are on opposite sides of the origin, so the stable 
invariant region includes the values £i < £ < 0. If the data can be set up in such 
a way that j/o and qo lie in the interval (£i,oo) for all xq, then the corresponding 
solution would have finite gradients for all time, and these would be nontrivial shock- 
free solutions. This topic is part of the authors' ongoing research. 

6. Generalized R/C structure. The R/C structure at a single contact discon- 
tinuity or entropy jump for the compressible Euler equation is fully analyzed in |14j . 
This is an analysis of how the R/C character changes when a wave crosses a jump 
discontinuity, where the system is a p-system on either side of the jump. Using this 
R/C structure on an entropy jump, a class of time-periodic solutions in compress- 
ible Euler equations with polytropic ideal gas has been studied in a series of papers 
[T4l [T5l [16] . In this section, we show that the R/C characters we have defined for 
smooth solutions are consistent with those results. Furthermore, by the study of R/C 
structures in the generalized Euler equations (|1.2p . (|1.3p with smooth and piecewise 
smooth entropy fields, we expect eventually to see a large class of time-periodic or 
quasi-periodic solutions as in |14[ 115) with both piecewise smooth and smooth entropy 
profiles. 
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We first consider how the R/C character can change at a fixed point in a smooth 
solution. 

Lemma 2. If {^-) v < (or equivalently {^f)h < 0j ; the backward R/C character 
can only change from C to R ( resp. R to C ), if the solution is forward C ( resp. R ); the 
forward R/C character can only change from C to R (resp. C to R), if the solution 
is backward R (resp. C). If > (or (^r)h > 0)> all the above R/C character 

changes will only happen in the opposite direction. 

Proof. Suppose the backward R/C character changes at (xo,io), then /3 = at 
(x 07 t ). Then by (|1.9[) . we have 

2 c 

which has the sign of a (^jr)h- Thus, the backward wave changes from C (j3 < 0) to R 
(f3 > 0) at (x ,t a ), if and only if a and {^f)h have the same sign. Thus if {^f)h < 0, 
then a < 0, and the crossing forward wave is necessarily compressive by Corollary [TJ 
All other cases follow similarly, and using Corollary [5] for the equivalence of derivative 
conditions completes the proof. □ 

In order to compare the R/C structure for smooth solutions to that of a single 
entropy jump, we recall the relevant argument from |14j . Discontinuities in weak 
solutions are governed by the Rankine-Hugoniot jump conditions, which are 

cr[v] = —[u], and a[u] = \p], (6-1) 

plus a third equation for energy conservation. Here a is the speed of the discontinuity, 
and brackets denote the jump [/] = fa — in / across the discontinuity. An entropy 
jump has zero speed, a = 0, so that (|6 . 1 1) reduce to [u] = = [p], that is 

u R = u L , and Pr=Pl- 

Across the entropy jump, the R/C structure changes as follows, see [T4] : 

Lemma 3. For ^ < 1, the backward R/C character can only change from R to 
C (resp. C to R), if the solution is forward R(resp. C); the forward R/C character 
can only change from C to R (resp. R to C), if the solution is backward R(resp. C). 
For I s > 1, all the above R/C character changes will only happen in the opposite 
direction. 

Comparing these Lemmas, we see that the condition < 1 for an entropy jump 
should be consistent with the condition (^ E ) v < for smooth solutions. Recalling 
that the change in pressure at an entropy jump is zero, we can assume p x = 0, and 
treat the condition ^ < 1 as c x < 0. Thus it suffices to show that, if p x = 0, then 



(— )„ < & c x < 0. (6.2) 

Pv 



Since f x = f v v x + f x , if p x = 0, we have 



so that also 



Px 
Pv ' 



. Px . C v Px + C x p v 

CvV x \ Cjr — C v ~r Cjr — 

Pv Pv 
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On the other hand, 



(Px\ PxvPv PxPvv 2c(Cxpt; Px c v) 



"Pv'" (Pv) 2 {Pv) 2 

where we recall p v = —c 2 . Comparing these, it follows that 

,Px_s _2c _ 2 

Pv Pv C 

which proves (|6.2I) . Thus our continuous R/C character is consistent with that of a 
single entropy jump. 

7. Transverse flow of ID MHD. In this section, we consider the motion of 
a compressible fluid coupled a magnetic field H = (Hi, H 2 , H 3 ), which satisfy the 
Magnetohydrodynamic (MHD) equations in Lagrangian coordinates. In one space 
dimension, the MHD equations can be written in the non-conservative form 

dv du\ 
dt dx 

dH 2 dux du 2 

-^T+P H 2-^. P H i^~ 

dt ox ox 

dH 3 dux du 3 
pHz-r. pH x - 



dt dx dx 

dui d , 1 

-W + d-x {P+ 2 



(P+-MHI+H 2 )) 



0. 


(7.1) 


0. 


(7.2) 


0. 


(7.3) 


0. 


(7.4) 


0. 


(7.5) 


0. 


(7.6) 


0. 


(7.7) 



— - ^ u ff 1 

dt dx 
du 3 _ dH 3 

dt dx 
dS_ 

~dl 

where p is the density, v = p^ 1 is the specific volume, (u\, u 2 , u 3 ) is the velocity field, 
S is the entropy, and p = p(v, S) is the pressure. The permeability pig is a positive 
constant, and Hi is constant in the one-dimesional model, c.f. [9l 112). 

The mathematical structures of (|7.ip - (|7.7p are totally different when Hi is zero or 
a nonzero constant. Here we briefly consider the easier case Hi = 0, which physically 
means that the magnetic field is transverse to the direction of motion. In this case, 
(fTTT]) — (f7T4|) and (|7J| form a closed system [9]: 

dv dui 

= 0, 

= 0, (7.8) 

= 0, 



dv 


dui 


~dt ~ 


dx 


dui 


dp 


~~dt 


dx 




dH 2 




dt 




dH 3 




dt 




dS 
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where 

1 H 2 + H 2 

H 2 = vH 2 , H 3 = vH 3 , p = _p+-Mo - 9 — -• 

2 

Thus in the special case of transverse flow, the one-dimensional MHD equations fit 
into the framework of (|1.2I) . (|1.3|) . and our previous results apply directly. 

We analyze this further for a polytropic ideal (7-law) gas, so that the pressure is 
given by 

p = Ke s/c " iT 7 , 

with adiabatic gas constant 7 > 1, and K, c v are positive constants. Then we get 



1 m + H'i 

Axv^' 1 + A 2 v~ 2 , 



p = Ke s / C " w-t + -Mo 2 2 3 
2 



where we have set 



Ai(x) = Ke s{x)/C " and A 2 (x) = ^/xo(fl|(a;) + ff 3 2 (z)). 
In order to further simplify our calculation, we now assume 7 = 2, and set 
B(x) = Ai(x) + A 2 (x), so that p(v, x) = B(x) v~ 2 . 
We then calculate 



c(v, x) = ^pT, = V2^B(xjv~V 2 , 

and, from (| 1 . T[) . 

h(v,x)= cdv = 2V2 t/b{x)v- 1/2 . 

J V 

Now solving for v and recalling fi — x, we obtain 



and so also 



Note that 



/i 4 ft 3 
PiKv) = ^7-57 \ and c i h ^) 



64B(ji) y 16 B(fj) ' 



dP -{ h iP) = c(h,n) and ^-(h,fi) = 



dh ' dh c(h, n) 

as intended by the choice of h, see [17] • Next, we calculate 

Vjj_ _ _h B(n) 
c 4%)' 
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where B = so that 



a — u x + h x 



h B(x) 



and j3 = u x 



h B(x) 



4 4 E(x) 

here a and /3 are evaluated at the point (x,t), and fi — x. From (|1 . 12|) . we calculate 

-1 B(p) 



which also leads to 



and thus, by ([T30]) . flTT]) . 



7(/i,/i) = 



ElL-l 



y : 



4/B 



h B 

h T 

5 B 



80 £(/z) 3 / 2 

-1 5(<i) 
20 B(fx) 3 / 2 

and = 



h 5 ' 2 , 



h 3 ' 2 
4VB 



K + - — 



ft s 
5 B 



Finally, using (|1.15p ~ (|1.17[) . we calculate the coefficients 

3 h 1 ' 2 , IB,, 
H2= 8 7f' and ai = 40Bl h = 

where we have used ()7.9|) . and, after simplification, 



B 



2 8 5 



1 1 / B 
W5B KW 2 

B 6 B 2 



B 2 



1 



-2^=1, 



1 



£5/2 5 #7/2 



B 7/2 V2 11 2 10 5 2 n 5 2 



-6 



G(/z)/ 2 



where we have set 



V ; 6 B(fj,) 2 

wherever B(p) 7^ 0. 

We use these in (|5.ip to calculate the quadratic 



h 1 ' 2 



6GI 2 ±2I£+ 



(7.9) 



and the corresponding dynamical system can be analyzed as above. It is clear that 
the growth of the quantities y and/or q depends critically on their size relative to /. 

8. Euler flow in a variable area duct. Finally, we consider the compress- 
ible Euler flow through a duct of varying cross section a(x). In (spatial) Eulerian 
coordinates (i,t), this flow satisfies: 



a{ = 0, 
(ap) i + (apu) x = 0, 
(apu)i + {apu 2 ) x + ap x = 0, 
(apE)i + (apEu + apu) x = 0, 



(8.1) 
(8.2) 
(8.3) 
(8.4) 
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where E is the energy, p is the density, u is the velocity, p is the pressure and a = a(x) 
is the diameter of the duct at position x, [3J HI [TT]. For smooth solutions, the energy 
equation (|8.4I) is equivalent to 



Si + uS x = 0, (8.5) 

where S is the entropy [3j. By introducing Lagrangian coordinates (x, t), which satisfy 

dx = apdx — apu dt, (8-6) 
dt = dt, (8.7) 

it is easy to check the smooth solutions of (|8. 1[) — 1|8.3[) and (|8.5[) satisfy 

v t -u x = 0, (8.8) 
u t + ap x = 0, (8.9) 
St=0, (8.10) 

where v is the specific volume per unit cross-sectional area, 

„ _ 1 
ap 

By by (|8.2p . the right hand side of (|8.6[) is an exact differential, so there is no difficulty 
defining x. By (|8.6|) . (|8.7p . a(x,t) satisfies the identities 

<Xt = ua, a x = v a, (a)t=u'd, and (d) x —vd, (8-11) 

where we denote 

da(x) .. d 2 a(x) 

a = — — — , and a = — —rr— , 

dx dx 

these describing the changing shape of the duct in spatial coordinates. 
We again assume we have an ideal poly tropic gas, so that 

p = K e s/c " v" 7 =Ke s/c " (at))" 7 . 

with adiabatic gas constant 7 > 1, c.f. [3J. The (Lagrangian) sound speed is 



T-l . 7 + 1 S 



c = y ' —apt, = y ^7 a 2 ?) 2 e 
Following [M] , we make the change of variables 

2V7?7 . 7-1 



s ^ and z=/ ^ d0 =fViiIfl-^. (8 .i2) 
Jv a ^ 2 m 7 1 



It follows that 

2 



_ n 2-, 

p = K p a 7 m 27-1, anc [ 

c = c(z, to) = A c a 2 mz7-i. 



where i"T„, if p and i"T c are appropriate positive constants, see |14j . 
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For C 1 solutions, the Lagrangian equations (|8.8[) - (|8.10[) are equivalent to 

z t + u x = 0, (8.13) 

a 2 to 

T-i ap 

ut + mca 2 z x + 2 — m x —"fpa x = 1 (8-14) 

TO 

mt = 0. (8.15) 

We define 

t-i 7 — 1 7-1 
a = ii T + a 2 m-z^ -| a 2 m^z and 

7 

n _2=i 7 - 1 _i=i 

p ~ u x — a 2 mzj a 2 rn^z, 

7 

which are the gradient variables used for the poly tropic ideal Euler flow in [Tj, adjusted 
to account for a(x). By similar calculations as in the proof of Theorem [T] or in [T], we 
obtain 

a' = fci(fc 2 (3a + (3) + (a(3 - a 2 )) + fc 3 (a - fi) + A(x, i), (8.16) 
f3' = -^(^(3(3 + a) + (a/3 - /3 2 )) + k 3 {/3 - a) - A(x,t), (8.17) 

where 

A(x,t) = ^^^[oa-Ta 8 ] + ——————— TTiTn x z 2 cl~^ ci : 

y ' J 8K C 1 ' J 2 7 

and the coefficients are 

fcl = 7 — — r-R'r^"' -1 , 

2(7-1) 

, 7 - 1 -^=i , 

«2 = -m x za 2 , and 

77+I 

, 3(7 -l) 3 _z±i. 7-I _ x . 
K3 = mza 2 a — ua a. 



By (I8T31) - (I8~15|) . 



so that 



2+1 7 1 y — 1 

z — z t + cz x — -KcZ-i- 1 (f3 H m x za ~), 

7 



= - T-^m.za-^. (8.18) 

Xcz— 7 

We plug (|8.18l) into (|8.16[) . move terms including z' to the left hand side, and 

T + l 

multiply by z 2 ^ -1 ' on both sides, to obtain 

(z5F=T7a<) + —a 5-z 2 <t-i) 2 ' _ — fc 3 z 2 <t-i) 2' 
27 K c 

7 — 1 ^ 7 — 1 o T+l 

= fci(fc2(3a — ^-m x za 2 ) + (-^-m x zo "a — a ))z 2 <t-i) 
+ k 3 (a + ~^m x za~^)z' : &^) + A(x,t)z^ lT ) . (8.19) 
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Furthermore, by (|8.13p - (l8.15l) . we have 

u' = u t + cu x = cf3 + jpa x 



y — 1 *y — 1 y — 1 

= —ma 3 z czm, x a ~+jpa x . (8.20) 

7 

By (|8.20p and following the proof of Theorem [5J (]8.19[) can again be transformed into 

y' = i>+{y) =a + a x y - a 2 y 2 , 

where y is a gradient variable, and ao, a\ are functions of m, m Xl m xx , z, a, a, a, and 
u, and 

7 + 1 7+1 I 

a 2 = ^7 -K c z^-v > 0. 

2(7-1) 

Similar considerations for backward waves yield a similar equation for the gradient 
variable q, namely 

q~ = rp-{q) = a - ai q - a 2 q 2 . 

A blowup result similar to Theorem [3] can again be obtained for this system. We omit 
the details since the calculations are tedious but similar to our previous calculation, 
and ao and a\ are more complicated expressions. Note in particular that they depend 
on u as well as z, which is a new feature of this case. 
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